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Abstract 

We present a numerical study of the shape taken by a spherical elastic surface when the volume it 
encloses is decreased. For the range of 2D parameters where such a surface may model a thin shell of 
an isotropic elastic material, the mode of deformation that develops a single depression is investigated 
in detail. It occurs via buckling from sphere toward an axisymmetric dimple, followed by a second 
buckling where the depression loses its axisymmetry through folding along portions of meridians. 
For the thinnest shells, a direct transition from the spherical conformation to the folded one can be 
observed. We could exhibit unifying master curves for the relative volume variation at which first 
and second buckling occur, and clarify the role of the Poisson's ratio. In the folded conformation, 
the number of folds and inner pressure are investigated, allowing us to infer shell features from mere 
observation and/or knowledge of external constraints. 

1 Introduction 

Let's consider a thin shell of an elastic isotropic material, such as a beach ball, and deflate it. What 
would be its shape ? 

This question is not restricted to garrulous familial shores : fundamental and applied physics nowadays 
presents legions of easily deformable soft objects, and knowing what governs their shapes gives the 
powerful possibility of inferring mechanical properties from simple observations, without contact. Among 
these deformable objects, an increasing number derives from spherical symmetry, that is omnipresent at 
scales where surface effects overcome volume forces such as gravity. The numerical study presented in this 
paper discusses the shapes taken by spherical thin shells of isotropic materials when their inner volume 
is decreased by a significant amount. Such a systematic and quantitative study will help deciphering 
conformations observed in e.g. Soft Matter (lock-and-key colloids[T], multiwall capsules[2], particles 
design through evaporation[31 13]), galenics (encapsulation [5]), microfluidics (microtanks [^) or medicine 
(ultrasound contrast agents [J), under the action of an external pressure or other possibly isotropic fields 
such as concentration in evaporation/dissolution phenomena, which also shapes objects in Nature |5]. 

When an elastic spherical shell has its inner volume lowered, it first deforms through in-plane compression 
that respects the spherical symmetry. Then it undergoes a symmetry breaking in order to relax a high 
stretch energy into much lower bending energy, by reversion of a spherical cap (creation of an axisymmetric 
depression, or "dimple"). The onset of this sudden transition, or buckling, under external pressure was 
studied long ago by Pogorelov and LandauflUl ITT] . They showed that the dimple should nucleate over a 
critical outside/inside pressure difference scaling like Y^o (-^) (in what follows, we will refer to this latter 
quantity by ^PLandau), where d is the shell thickness, R its radius and Yan the Young modulus of the 
material that makes it up. Its edge (or "rim") has a transversal extension \/dR. One of the key of their 
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calculation being the assumption that buckling occurs for dimples such that maximum deflection is of 
order d, spherical geometry imposes then that y/dR is also the radius of the dimples that forms. Besides, 
classical buckling analysis provided dependence of the buckling pressure with the Poisson's ratio[51 [TUllllj . 

Results focusing on deformations through further deflation are, mainly, more recent. A few months ago, 
stability analysis allowed a detailed study of buckling toward axisymmetric conformations [T2| . Experi- 
mental [m [3l [m [iSl [m [H [TT] and numerical O [TH [l7] deflation studies showed shapes holding several 
dimples, also called "multiple indentation". These conformations compete with experimental observa- 
tions of shapes holding a single depression |14[ IT5] . possibly losing axisymmetryj^Ol or exhibiting folding 
perpendicularly to the rim |21[ [15] . Similar shapes are observed in shells under a point load [22l [23] 
or pressed against a wall [241 . Secondary buckling by folding of the single depression, also called 
"polygonal indentation", under isotropic constraint was numerically retrieved with surface models pnil8| . 
Thin shells with a single depression, either axisymmetric or polygonal, appear to present a conformation 
of lower energy than the metastable multiple indentation[5Tl [THl [T7]. In the case of an axisymmetric 
dimple, this can be easily understood since elastic energy mainly concentrates in dimple edges as bending 
energy, with an energy per edge length that weakly varies with dimple size. Hence dimples coalescence 
lowers the total elastic energy |25]. Nevertheless, more than one dimple may nucleate if the deflation is 
rapid enough, leading to metastable[3l [HI [17] multi-indented shapes. The term "rapid" is to be taken 
on a wide acception here. Experimentally, it may correspond to situations where dissipation (due to 
material viscosity or to fluid flows accompanying the deformation) prevents dimple growth, which favors 
secondary nucleation once l^Phandau is reached, and where subsequent kinetics prevents thermally acti- 
vated coalescence between adjacent dimples. Numerically, minimization may reproduce such metastable 
situations [3l[IHl[T7], since (i) large volume increments favors the creation of extra dimples, by making 
difficult to find the cooperative displacement of vertices that corresponds to rim rolling in dimple growth 
(ii) depending on the way curvatures are calculated, energy barriers that prevents from dimple coales- 
cence may be overcome or not. For "slow" deflations, a single dimple can appear and grow, or freshly 
nucleated dimples may coalesce into a single one. Such "slow" deflation provides an axisymmetric bowl- 
like shape, that may undergo under further deflation a transition toward a non-axisymmetric depression, 
i.e. polvgonal indentation|2H 118]. 

We present here a systematic numerical study of such "slow" deflations leading to shapes with a single 
depression. In this purpose, we used a surface model taking into account recent developments, presented 
in section [2l We clearly expose the correspondence between 2D parameters of the model surface, and 
3D properties of the real object of nonzero thickness, expliciting the role of the different significative 
parameters. Particular emphasis is put on a parameter often underconsidered: the Poisson's ratio. 

The whole study allows to determine parameters of importance for the transitions sphere — > axisymetric 
bowl (section [3]) and axisymmetric bowl — > polygonal indentation (section [4]), both for the detailed 
shape in polygonal indentation, and for inner pressure. Furthermore, we took particular care to provide 
empirical dependence laws for practical use. 

2 Surface model 

Surface model, where out-of-plane and in-plane deformations are formally uncoupled, is for long consid- 
ered as valid to describe the deformation of thin sheets (plates or shells) [11] [5^ . For thin sheets without 
spontaneous curvature {i.e. an elementary surface portion of the sheet, freed from constraints exerted by 
surrouding material, remains flat at equilibrium) , the energies per surface unit that are to be considered 
in this surface model are of two kinds: firstly, a curvature term that can express ^kc^ -I- 'Kg |271 111) , 
where c = -I- and 9 — ^ ^ ^ ^"^^ respectively the mean and Gaussian curvatures (i?i and R2 
being the local principal curvature radii) , and k and k are respectively the mean and Gaussian curvature 
constants [28). The other term may be written, in a Hookean linear model: ^^ijKijkieki, where Sij and 
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Kijki respectively represent the two-dimensional strain and elasticity tensors for in-plane deformations. 
For an homogeneous and isotropic surface, the nonzero terms of the two-dimensional elasticity tensor 

Kxxyy — -^yyxx — ~^~i7r^ and Kxyxy — ^yxyx — TTT^TT : with Y2D the 
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^yyyy — i-u^^ ' ^^^^yy ~ '-^yyxx — i-i,'^^ ^'-xyxy — j-^yxyx — i+y^^- 

two-dimensional Young modulus and V2d the two-dimensional Poisson ratio, which is comprized between 



-1 and 1 This in- plane elasticity term can be rewritten as 
of concision. 
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In a linear approximation, the relation between the 2D parameters and the 3D features of the plate 
(Young modulus Ysd, Poisson's ratio v^d, thickness d) with zero boundary tangential constraints [11] is 
expressed as (detailed e.g. in PH or [7]): 



V2D = VzD = V 

YiD = YzDd 
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Y3D 



12(1 + 1^) 



(1) 
(2) 
(3) 

(4) 



Since for bulk materials the maximum value of 1^30 is ^ for thermodynamic reasons [TT], one can notice 
that the range of 2D Poisson's ratio that effectively describes a thin plate of an isotropic material is limited 
to a maximum value of i. In other terms, even a thin plate of an incompressible isotropic material cannot 
behave as an incompressible surface (where 1^20 = 1), thanks to the possibility of having its thickness 
varied. On the other limit, Poisson's ratio can reach -1 as a lower value, but negative values correspond 
to less common "auxetic" materials. 

Conversely the thickness of the plate, as a function of 2D parameters, writes: 



d= JI2 (1-1.2) 



(5) 



For describing surfaces with asymmetric properties, the notion of "spontaneous curvature" was introduced 
by W. Helfrich It was recently shown that to describe the deformations of an initially stress- free 

thin shell of radius R, the three contributions (in-plane, mean curvature and gaussian energy, that may 
have a non- vanishing part even for closed surfaces, depending on the definition of a nonzero spontaneous 
curvature) can be rewritten in an easily computable way as [29) : 
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with Cq = being the effective spontaneous curvature, and f — ^"2/?^'''" effective surface tension. 

This expression is slightly different from the one used in [21J, hence we will quantitatively discuss, in the 
results, modifications induced by the use of this more complete expression. 

More generally, we will consider the influence of sphere size through the use of the adimensionalized 
Foppl-von Karman number [50] 7 = , that gives the order of magnitude of the ratio between in- 

plane and out-of plane deformation energies An elastic surface with the energy given in equation 
([6]) can effectively describe a thin shell of an isotropic material if 12 (l — u"^) /7 ^ 1, in addition to the 

condition 1^20 < 5- In this range, 7 roughly scales like (-f )^. Out of this range, such a surface model does 
not correspond to any thin shell of an isotropic material ; it can nevertheless describe different types of 
objects, e.g. gel-phase vesicles [T7| that can hence be considered as thin shells of non isotropic materials. 
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Numerical experiments are performed by minimizing the elastic energy as expressed in equation ([6]) for 
different inner volumes, with the free software Surface Evolver[5T]. A whole in-silico deflation experiment 
(from Vinit = l""^^ to w 0.2 x Vinit) is realized through a succession of different equilibrium states, 
these latter found according the process described in [21], and calculated successively for inner volumes 
decreased by steps of at maximum 2% of the initial volume (steps amplitude is reduced in some situations 
in order to avoid the nucleation of secondary dimples). 

3 First-order transition toward axisymmetric depression 

Deflation of a spherical elastic surface, at imposed either volume or external pressure, causes an abrupt 
buckling from the spherical conformation in order to release in-plane compressive stress. The purpose 
here is to compare the numerical approach described in the previous paragraph to known features of this 
buckling, for the range of parameters that scans the generality of thin shells of isotropic material. 

3.1 Buckling pressure 

First buckling relaxes in-plane constraints, and causes drastic drop of the inside/outside pressure differ- 
ence AP = Pext — Pint- Figure [T] displays typical evolutions for the pressure difference: first a linear 
increase followed by a drop at first buckling, after which pressure difference varies in a much lesser ex- 
tent. Linear behaviour is expected before the first buckling due to the relation between pressure and 
elastic energy P^xt — Pint = dEeiastic/d (AF) (detailed in [7|), and quadratic dependence of Eeiastic with 
= V,mt " V (see e.g. HH). 

In the simulations presented here, designed not to be stuck in multi-indentation conformations of higher 
energy, or other less stable, the first buckling leads to a single axisymmetric dimple. It is expected to 
happen when the external overpressure AP = P^^t — Pint reaches the critical value [HI [T^] : 

AP. = 2[3(l-.^)]-^/^xr3z,(|)' (7) 

As expected, first buckling in our simulations effectively occurs at a pressure difference of order APj^andau = 
YsD (■^) (section[T]). For a given Poisson's ratio, the incertitude due to discrete volume increments does 
not allow to conclude that APbuckHngi / AP^andau is affected by 7 (see e.g. figure [T] displaying several 
cases at ly = —0.5 ). The effect of the Poisson's ratio is displayed on figure [5] it shows that the pres- 
sure which induces buckling in our simulations quantitatively follows the theoretical equation [71 which 
reinforces the validity of our approach. 

3.2 Buckling volume 

The mechanism of buckling in axisymmetric conformations was quite recently investigated in detail by 
Knoche et aZ|12|. with the study of various metastability branches. At first significative order, their 
calculations show that for the trivial isotropic ("spherical") deformation, the relative volume variation 
due to an external overpressure A/^ — Pext — Pint expresses: 

AV"_3(l-;^) R AP 
~V~ ~ 2 7 1^ 

Hence the deflation at buckling pressure writes}^ I12j: 



4 




Figure 1: Pressure difference AP = Pext — Pint adimensionalized by APLandau = Ysd {^) (or 

1 /2 

[12Y2DK (1 - 1^'^)] /R^ in 2D parameters), as a function of the relative volume variation, for surfaces 
of similar Poisson's ratio v = -0.5. Dotted green: 7 = 1.17 x 10^ = 8.8 x 10^^); blue: 7 = 6.07 x 10^ 
(-| = 1.22 X 10~2), interrupted orange: 7 = 3.22 x 10^* (|r = 1.67 x lO'^)^ magenta: 7 1.17 x 10^ 
(| = 2.31x10-2), light blue: 7 = 9.33 xlO^ (| = 3.11x10-2), parma: 7 = 4.67x10^ (| = 4.39 xlO^^). 
Notice scale switch at = 0.2. Points (corresponding each to a minimization) emphasizes two curves 
with typical behaviour after buckling: increasing (green curve), and plateauing (orange) . 
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Figure 2: Pressure i^P)i,ucki ~ (Pf^-^t ^ -Pi^-t) buck i &i^st sphere— >-bowl buckling, adimensionized by 
^Phandau- Error bars represent the range of this quantity for 7 between 4.67 x 10"^ and 2.33 x 10^. Line: 
theoretical value 2 [3 (l — t^^)] from equation[71 
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Figure [3] presents values of (■^t") from numerical simulations, as a function of a combination of 
7 and u translated in 3D parameters: here also the theoretical equation [5] is quantitatively retrieved. 
One may notice (since by essence {^-)^,^^f.l ^ 1) that the spheres-bowl transition vanishes for > 
^y{l + ly)/3(Y—u) , which indicates a destabilization of the axisymmetric bowl for the most auxetic 
materials. This was qualitatively expected from the k divergence in this limit, which makes curvature 
deformations prohibitive compared to in-plane compressions [32] . 

Using the 2D parameters of the surface model, relation |S] also writes: 

= 6(l-j/)7~^ (9) 



AV 



buck 1 



Extending our purpose out of the range of 2D parameters that effectively can describe a thin shell of 
isotropic material, we may remark that extrapolation to v = 1 induces vanishing of {^~) i^y^^i^ hence 
destabilization of the spherical conformation : this limit corresponds to incompressible surfaces, that will 
necessarily undergo a deformation implying curvature even for the smallest volume decreases, since area 
variation is prohibited. For the thickest shells, i.e. 7 < 80, shape transition occurs not any more through 
sudden inversion of a spherical cap, but by slowly deforming into an ovoid, that flattens at the location 
of future depression under further deflation. We did not, here, specifically study this extreme behaviour. 

For practical purposes, we may notice that the large range of parameters explored shows that in terms 
of volume, the onset of buckling mainly depends on the relative thickness 4, with only a weak influence 
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Figure 3: Relative volume variation {^')f,y^^f^i ^-t the sphere— >-bowl buckling (lower and upper limit of 
the error bars are indicated respectively by upward and downmard triangles). White: v — 0.5 ; light 
gray: i/ = ; dark grey: ly = —0.5 ; black: i/ = —0.8. Grey line: = 6 (1 — v)^^^ (this is 

equation |S] expressed in 2D parameters). Dotted line: relation established in for v = 1/3: taking 
into account the gaussian curvature elastic energy increases volume threshold by at maximum 20% for 
the thinnest shells of the shells previously studied. Illustrative inserts display deflated spherical surfaces 
at 7 = 4666 and v = —0.8 (hence = 0.0304), with relative volume variations respectively = 0.161 
(spherical) and 0.167 (buckled, section view, same scale). 
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Figure 4: Left: schematisatioii of a spherical surface near complete deflation. Inner part of the rim 
endures a compressive stress in the direction indicated by the red double arrow. Relaxation of in-plane 
deformation occurs, for the thinnest shells, via undulation deformation, generating wrinkles (folds) all 
along the rim (examples schematized with blue interrupted lines). Length I stands for lateral extension 
of the wrinkles. Right: simulation with 7 = 9.33 x 10^, v = 0.5 (-| = 0.031) and ^ 0.562 ). Maxima 
of undulation are stressed with blue interrupted lines ; circle locates the zone of high curvature that forms 
the apex of the s-cone. 

of the Poisson's ratio. Since this latter ranges between and \ for common materials (i.e. non auxetic, 
with > 0), the prefactor of varies between \fZ and 1, which is much narrower than the range in 
that can be explored. 

4 Second- order transition toward polygonal depression 
4.1 Location of the transition 

In the axisymmetric bowl shape, global bending of the rim on the equator costs in-plane deformation: 
extension on the outer side of the rim, and compression on the inner one. For the thinnest shells, compres- 
sive stress parallel to the equator leads to a secondary buckling, where the inner side of the rim undulates 
to adapt to axial compression (fig. 31 left), forming folds, or "wrinkles", that deform the axisymmetric 
depression into a roughly polygonal shape (fig. |4l right). Such a conformation mainly involves curvature 
deformations, much less energetic |24| than compression energy that quadratically increases with 
Fig. [5] shows how elastic energy dispatches between in-plane and out-of-plane deformation energies in 
a typical numerical deflation. Wrinkles match with the rim through a zone of high curvature that has 
a folding role similar to what realizes the apex of d-cones |331 [Ml [55] , except that the surface is not 
developable but is spherical, hence the concept of " s-cones" proposed by Reis and Lazarus [36] . 

Secondary buckling from axisymmetric bowl shape to polygonal indentation, quite smooth, is harder to 
detect than the first one (fig. [5]). The corresponding relative volume variation (■^t")^^;,^. 2 determined 
on one hand by the maximum deflation before loss of axisymmetry, and on the other hand by the at 
which the rim presents convex zones under axial observation (as shown on lower part of Fig. [S] subfigure 
c). Figure [5] presents a typical shape phase diagram for v = —0.5, with three distinct zones: spherical 
coformation, axisymmetric bowl and wrinkled depression. One may notice that for the thinnest shells, 
the incremented deflation we numerically performed shows a direct transition from the sphere to the 
wrinkled bowl. Relative volume variation at first buckling nevertheless obeys a single power-law on the 
whole range of relative thicknesses. We focus in this section on the second buckling, hence considering only 
deflations where wrinkles appear on an already existing axisymmetric conformation. Figure [7] displays 
how (■^7~)(,„cA;2 ^^''i^^ with the Foppl-von Karman number 7 for different Poisson's ratios. Data indicate 
a dependence in 7, but scattering prevents from concluding on an influence by v ; linear regression on 
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Figure 5: Deformation energies at different relative volume variations (volume step « 0.6% of the initial 
volume), adimensionalized by the curvature constant k, for 7 = 9.33 x 10'^ and v ~ 0.5 (-^ = 0.031). 
Interrupted line: Hookean in-plane deformation energy, linked to the second significative term of eq. ([6]), 
i.e. ^€ijKijki€ki. Grey continuous line: elastic energy linked to the first significative term : (c — Cg) . 
This term does not vanish at = because Cq ^ Cq ; in the total elastic energy, it is counterbalanced 
by the est and the effective surface term expressed by 7e// in equation Effective surface energy 
(not represented here) varies at maximum by — 0.33k , and on an amplitude 0.03 k in the nonspherical 
conformations. Grey zone indicates the second transition from axisymmetric bowl shape to depression 
with inner wrinkles ("polygonal indentation" [18]), determined as explained in the text. First transition 
from sphere to axisymmetric bowl occurs through abrupt decrease of the in-plane deformation energy, at 
^ = 0.03. 
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logarithms [37] provides with a correlation coefRcieiit of -0.99 
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= 8470 X ^-^■'J^'^ (10) 

buck 2 



This expression is compatible with previous results obtained for v = l/i with slightly different numerical 
models: the shell without gaussian curvature evoked in section r2|21| . and the spring model without 
spontaneous mean curvature of reference |18) (both are represented on figure [7] for their range of validity). 

Extrapolating equation [TO] up to = 1 suggests that this secondary buckling does not happen, i. e. 
single indentation keeps its axisymmetry, below a threshold value "fc,buck2 = 4170. In tridimensional 
parameters, relation (|10[) expresses: 



V ^ / buck 2 VVl - V^J 

The dependence in i/ for non-auxetic materials is even weaker than for the axisymmetric buckling since 
it plays at maximum by a factor 4/3. Besides, for the range of Poisson's ratio studied here, there is no 
axisymmetric conformation to be expected for d/R < 0.003. Relation ([TT]) also implies that wrinkles are 
not expected when d/R > 0.054 x y/l — : a particular consequence is that in wrinkles prevention, a 
very auxetic material (with — > — 1) may help. 



4.2 Characterization of buckled shapes with the number of wrinkles 

For the thinnest shells that undergo polygonal indentation, the most conspicuous feature is the number 
W of wrinkles, or s-cones. Figure [S] shows the evolution of W for a typical deflation: first does W decrease 
while the freshly nucleated and still very flat depression hollows and enlarges, then it increases again. 
Data are quite scattered: there is a typical noise of order ±1 on W, that has no observable correspondence 
in smooth energy curves (fig. [T])- In order to decrease data scattering, we calculated Wdefiated as the 
average value of W between = 0.53 and — 0.76 (these values have been choosen in order to 
cover, for all the simulations, a significant range of relative volume variation before autocontact, this 
latter happening around ~ 0.9). Values of Wdefiated are comparable with results from the previous 
model, which did not take gaussian curvature into account, and indicates a scaling law in 7^1/^ [211 [7]. 
The dependence of Wdefiated with 7 and v is shown on figure [^ expressed in 3D parameters. It shows 

a scaling in (-^) , which provides clues on the typical transversal size / of the s-cones (presented on 
figure [3] left). Since s-cones stand alongside one another on a length which is of the order of an equator, 

we can estimate I as wT^JTTl' ^^^'^^ best fit of Wdefiated with dependence in(-^) (figure [H]) can be 
expressed as: 

I ^ 6.7 VdR (12) 

This result is fully comparable to the wrinkles wavelength 4.7-\/ dR that can be calculated from recent 
results by Vella et al [TO] on the indentation of strongly pressurized shells. As shown in reference |11| . 
V dR arises naturally from balancing the bending and in-plane deformation energy of a small deformation 
on a spherical shell. Recent results by|12j showed that VdR scales also for the rim width even in large 
axisymmetric depressions ; we confirm here that it governs also other types of large deformations such 
as s-cones transversal size in polygonal depressions. 



5 Other postbuckling features 

First buckling and its consequences were presented in section [3] For application purpose, we may in- 
terest to the inside/outside pressure difference after its drop at first buckling. Careful examination of 
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Figure 7: Points: relative volume variation at which the axisymmetric depression becomes "polygonal", 
versus ^Ylj^ (which is to be identified with djR when = 0). White: v = 0.5 ; light grey: v — Q \ dark 
grey: v ~ —0.5 ; black: v ~ —0.8. Wide grey line: equation [TT] Continuous lines indicate the location 
of the first-order buckling for different Poisson's ratio (from equation [8|) , under which only the spherical 
conformation is to be found. Non-continuous lines: in their domain of validity, equations describing the 
axisymmetric/wrinkled bowl transition published in previous works for = 1/3 (see text): dotted black 
refers to [3T], interrupted grey to |18j . 
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Figure 8: Number of wrinkles (s-cones) held by the single depression after the secondary buckling of 
figure [5l Poisson's ratio v = 0.5; Foppl-von Karman number 7 = 6.06 x 10'*(-^ = 0.0122). Lower part: 
conformations at points indicated on the main figure. 
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Figure 9: Number Wdefiated of wrinkles (s-cones) at the end of a numerical deflation (averaged between 
= 0.53 and = 0.76). error bars are taken as the standard deviation on this range, with a minimum 
value of ±0.5. White: v = 0.5 ; hght grey: v = dark grey: v = —0.5 ; black: v = —0.8. Continuous 

line: Wdefiated = OMO x (|)"^ 



the numerical data revealed supplementary features, uncorrelated with the occurrence of the secondary 
buckling previously exposed : two types of behaviour clearly appear in the evolution of the reduced 
pressure ^p^^^ — during deflation. 

For ^ > 0.014, pressure difference AP presents the type of evolution calculated by [12], *.e. quasi- 
plateauing after buckling (variation of about 15% during the whole deflation, plus some occasional dis- 
persion due to numerical procedure), up to autocontact. Furthermore, an order relation is respected: at 
every volume step, the ratio ^J^^ ^ — weakly increases when 7 decreases. This is to be observed in fig. 
[1] for the 4 curves corresponding to the highest relative thicknesses. 

For the thinnest shells (-^r < 0.012 in the simulations performed), ^J^^ ^ — regularly re-increases with 
deflation after the pressure drop, crossing successively the curves at smaller 7's (as shown on the 2 
"thinnest" curves of flg. [ij . 

In order to extract a general behaviour from these differents observations, we focused on ^p^'"'^' ; the 
minimum value of ^p^^^ — after buckling. Figure [TOl shows that the cross-over between the two regimes 
around « 0.014 also corresponds, for each to the minimum of relative pressure drop after the first 
buckling. In the plateauing regime, _^p^""^" stands for the plateauing value for ^ > 0.014 ; it shows 

a power-law of the type ^p^"*'^' = a{v) x (■^)°^. Similarly to what was done in section[21 and since 
a {v) appears to be even (curves aX v = 0.5 and v = —0.5 almost mix up on figure [TU]) . we looked for a 
prefactor of the form (l — v^^ , minimizing x for the best fits at -^r > 0.014. This led us to propose the 

master curve presented on figure [Tl] of formula ^y^^" ^ ~ '^^)^ '^'^^ = 0-75 x (-^)^^. We do not have 
for the moment theoretical clues to justify these two successive fitting operations, but (i) it allows to 
describe numerical results in a very condensed way for > 0.014 , and (ii) for all the shells these reduced 
values impressively gather on a single curve, for 7 ranging from 8 x 10^ to 4.7 x 10^, and for v between 
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Figure 10: Minimum inside/outside pressure difference after buckling (cf figure [T]), adimensionalized by 
Landau pressure, logarithmic representation. Black: v = 0.5; dark grey: i/ = 0; light grey: v = —0.5; 
white: v = —0.8. Interrupted line separates the two types of evolution of the pressure after the first 
buckling: increasing or plateauing (see figur£[T]). Black line indicate slope 0.5. 



-0.8 and 0.5. This result, exposed on figure [TT] using 3D parameters, is expected to be of practical use 
for all experiments involving defiation controlled by the volume. On a more conceptual point of view, 
plot clearly confirms two different scalings of the pressure during deflation, around a threshold in relative 
thickness « 0.013. This may be an indication of the existence of different ways to accomodate 

s-cones on a sphere, and requires further investigations. 



6 Conclusion 

Systematic numerical study of the buckling of a spherical shell, in the conformation with a single depres- 
sion, allows to sketch the influence of the different geometrical or elastic parameters through quite simple 
theoretical or phenomenological laws. The surface model can be translated in 3D parameters, that are 
the shell's thickness, and the two elastic parameters of the material that compose it: Young modulus and 
Poisson's ratio. 

At imposed volume, the Young modulus does not play on the shape. Results showed that the first 
transition (toward axysimmetrically buckled shape), and the second one, with appearance of wrinkles, or 
"s-cones", is mainly driven by for non-auxetic (i.e. with positive Poisson's ratio) materials. For auxetic 
materials, the Poisson's ratio may have a determining importance, by strongly displacing transitions 
toward higher values of the relative volume variation, up to possible vanishing. Decreasing the Poisson's 
ratio down to very negative values stabilizes spherical defiation at the expense of dimples creation, and 
axisymmetric dimples against appearance of wrinkles. 

The number of wrinkles indicate a dependence in (■^) , that confirms ^/ dR as the accurate scaling 
for elastic deformations of elastic spherical surfaces. 
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Figure 11: Pressure master curve after first buckling: x (l — i^^) ' versus Black squares: 

V = 0.5; black diamonds: v — 0; light grey diamonds: v = —0.5; white diamonds: u = —0.8. Continuous 
line: ^x(l-.^)°-^^^ = 0.75 X (1)^-1 

The Young modulus scales pressure features: critical inside/outside pressure difference that triggers first 
buckling, and plateauing pressure after buckling. Detailed behaviour, that is shown to reduce to a master 
curve, opens the possibility for two different wrinkling regimes. 
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